January 20, 2012, Ver 0.8a 

Preprint typeset using WT^K. style cmulatcapj v. 5/2/11 



LOCAL ANISOTROPY, HIGHER ORDER STATISTICS, AND TURBULENCE SPECTRA 
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ABSTRACT 

Correlation anisotropy emerges dynamically in magnctohydrodynamics (MHD), producing stronger 
gradients across the large-scale mean magnetic field than along it. This occurs both globally and 
locally, and has significant implications in space and astrophysical plasmas, including particle scatter- 
ing and transport, and theories of turbulence. Properties of local correlation anisotropy are further 
documented here by showing through numerical experiments that the effect is intensified in more 
localized estimates of the mean field. The mathematical formulation of this property shows that local 
anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates 
to higher order correlations can be understood in connection with the stochastic coordinate system 
inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to 
higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus 
establish that the local structure function is not a measure of the energy spectrum. Evidently the local 
enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon 
must be understood in terms of higher order correlations, fourth-order and above. 



1. INTRODUCTION: CORRELATIONS AND ANISOTROPY 

Correlation and spectral anisotropy play important 
roles in current studies of solar wind and astrophysi- 
cal plasmas, having significant impact on descriptions 
of the turbulence cascade, particle scattering, the nature 
of kinetic dissipation, and the transport of turbulence. 
Evidence from experiments, numerical simulations, the- 
oretical arguments, and spacecraft observations has con- 
sistently supported the conclusion that magnctohydro- 
dynamic (MHD) turbulence leads to states characterized 
by gradients that are relatively stronger when measured 
perpendicular to the large-scale magnetic field and rel- 
atively weaker when measured parallel to it (Robinson 
& Rusbridge 1971; Zweben et al. 1979; Ren et al. 2011; 
Shebalin et al. 1983; Bondeson 1985; Carbone & Veltri 
1990; Oughton et al. 1994; Montgomery & Turner 1981; 
Matthaeus et al. 1990; Bieber et al. 1996). This corre- 
lation anisotropy has been quantified both globally and 
locally, by varying the definition of the mean magnetic 
field. The local form, being of greater magnitude, has re- 
ceived substantial attention in recent years (Cho & Vish- 
niac 2000; Milano et al. 2001; Horbury et al. 2008; Tessein 
et al. 2009; Podesta 2009; Luo & Wu 2010; Wicks et al. 
2010, 2011; Chen et al. 2011). The present paper will fo- 
cus on the relationship between these forms of anisotropy, 
providing a better characterization of the scale depen- 
dence of local anisotropy. A particular conclusion will 
be that correlation anisotropy affects statistics at all or- 
ders, including but not limited to the energy spectrum 
and other second-order statistics. We also present evi- 
dence that the enhancement of local anisotropy over the 
global value is due entirely to higher order statistical ef- 

1 Department of Physics and Astronomy, University of 
Delaware, Delaware, USA 

2 Dipartimcnto di Fisica, Univcrsita della Calabria, 1-87036 
Cosenza, Italy 

3 Departamento de Fisica (FCEN-UBA & IFIBA-CONICET), 
Buenos Aires, Argentina 

4 Department of Mathematics. University of Waikato, Hamil- 
ton 3240, New Zealand 



fects. 

The basic physics leading to this strong perpendicu- 
lar spectral transfer and enhancement of perpendicular 
anisotropy has been elucidated in the context of incom- 
pressible homogeneous MHD turbulence (Montgomery 
& Turner 1981; Shebalin et al. 1983; Bondeson 1985; 
Grappin 1986; Oughton et al. 1994). For this model, all 
modal couplings are triadic (involving wavevectors k%, 
k2, ks, such that k\ + &2 + &3 = 0), but in the pres- 
ence of a strong uniform DC magnetic field BqZ, only 
those couplings that have one (or all three) wavevec- 
tors perpendicular to z will proceed unattenuated by 
propagation effects. All other couplings are suppressed 
to a progressively greater degree as B is strengthened. 
The appearance of stronger perpendicular gradients im- 
plies that spectral transfer to the perpendicular wavevec- 
tors is more robust than that to the parallel wavenum- 
bers. This corresponds to appearance of characteris- 
tic spectral anisotropy, with enhancements of high k± 
power relative to k\\ power, an effect that become pro- 
gressively stronger at smaller scales (Shebalin et al. 1983; 
Oughton et al. 1994). This organization of the spectrum, 
is built into models such as Reduced MHD (RMHD) 
(Kadomtsev & Pogutse 1974; Strauss 1976; Montgomery 
1982) and related models such as the Goldreich & Srid- 
har (1995) steady-state (GS) model, and the nearly in- 
compressible MHD (NIMHD) model (Zank & Matthaeus 
1992). Indeed strong anisotropy of MHD correlations is 
observed in laboratory devices (Robinson & Rusbridge 
1971; Zweben et al. 1979; Ren et al. 2011), in the solar 
wind (Matthaeus et al. 1990; Tu & Marsch 1995; Bieber 
et al. 1996) and the corona (Armstrong ct al. 1990), and 
is an apparent requirement for consistency of scattering 
theory with solar energetic particle (SEP) observations 
(Bieber et al. 1994; Droge 2005). 

However anisotropy of the energy spectrum is not the 
only effect associated with correlation anisotropy and en- 
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hancemcnt of perpendicular gradients. 5 For example, 
one can discuss the shape and orientation of structures, 
such as "eddies" (Cho & Vishniac 2000) or current sheets 
(Dmitruk et al. 2004). It is well known that one of the 
powerful nonlinear effects of turbulence is the produc- 
tion of coherent structures that are progressively smaller 
at small scales. This leads to intermittency as mea- 
sured by nonGaussianity and higher order statistics, as 
seen in hydrodynamic and MHD models (She & Leveque 
1994; Politano & Pouquet 1995). Simulations show that 
such nonGaussian statistics are generated very rapidly by 
the cascade of excitations to smaller scales (Wan et al. 
2009). In MHD, under a wide variety of conditions, 
the associated coherent structures take the form of en- 
hancements of electric current density in sheets or fila- 
ments (Matthaeus & Lamkin 1986; Carbone et al. 1990; 
Veltri 1999). When the turbulence is anisotropic rel- 
ative to a large-scale field (e.g., Dmitruk et al. 2004), 
the current structures tend to align in that direction as 
well. This suggests that higher order statistical quanti- 
ties (fourth-order moments, kurtosis, etc) also must be- 
come anisotropic. 

The above reasoning leads to an expectation that 
statistics at all orders might be involved in correlation 
anisotropy, but many "theories of MHD turbulence" dis- 
cuss exclusively the properties of the second-order statis- 
tics; that is, the energy spectra and associated correlation 
and structure functions. This emphasis follows naturally 
from the use of wavevector spectra and two-point, single 
time correlation functions as measures of the distribution 
of energy in spatial structures of varying size (Monin Sz 
Yaglom 1971, 1975), as exemplified by the classic Kol- 
mogorov (1941) theory. Recognizing this background, 
it is not entirely surprising that recent interest in local 
anisotropy has sometimes focused on an interpretation 
of this effect as a local energy spectrum (Cho & Vishniac 
2000; Chen et al. 2011; Wicks et al. 2011). The discussion 
below provides in effect a critique of this interpretation. 

To render the discussion specific, we assume a statis- 
tical description that is homogeneous in space and sta- 
tionary in time. The brackets (• • • ) denote an ensemble 
average, which by invoking a classical ergodic theorem 
is approximated in practice by a space or time average. 
It may be possible to define the statistical ensemble in 
other ways, but here a classical statistical framework is 
assumed (Monin & Yaglom 1971, 1975). The two-point 
correlation function measuring the statistical relation- 
ship between the magnetic field fluctuation at points x 
and displaced position x' = x + r is also the Fourier 
transform of the wavevector energy spectrum, 6 S a p(k). 
The definitions are (e.g., Monin & Yaglom 1971, 1975) 

R a p(r) = (b a b' p ) = J d 3 k S aP (k)e lk - r , (1) 

where we abbreviate b a = b a {x, t) and b'p = bp(x + r, t), 

5 Here we are concerned with correlation (or spectral) anisotropy, 
e.g., anisotropy of |b(k)| for varying k. Anisotropy of the variance, 
i.e., direction of b(x), also occurs, see Belcher & Davis (1971), but 
is a distinct issue, which however may become linked in certain the- 
ories, e.g., Zank & Matthaeus (1992); Goldreich & Sridhar (1995). 

6 For three dimensional isotropic Kolmogorov turbulence the 
omnidirectional energy spectrum is £{k) = 4Trk 2 S aa (k) (sum on 
repeated indices), with £(k) ~ A; -5 / 3 in the inertial range. 



and suppress the time argument. 7 Another two-point 
statistic is the second-order structure function 

D(r) = (\b-b'\ 2 ) =25b 2 -2R(r), (2) 

which is obviously related to R a p . Here Sb 2 = R aa (0) is 
the total variance of the fluctuations and R(r) = i? Qa (r). 

Existence of a homogeneous statistical ensemble im- 
plies that there exists a probability density V(T) that de- 
scribes fully all realizations of the turbulence, where each 
realization is labeled by an index T. This density may 
be projected onto the space of two-point statistics, arriv- 
ing at a probability distribution function (pdf) P 2 (b,b') 
that characterizes all of the two-point statistics of any 
order. Thus one may equivalently express the structure 
function in two ways as 

D(r)=J^dT [\b-b'\ 2 ] r V(T) (3) 
= J J d 3 bd 3 b' \b-b'\ 2 P 2 (b,b'), (4) 

where — 6'| 2 ] r is the value of the square- bracketed 
quantity for realization T. Similar relations hold for the 
spectra and two-point correlation functions. Note that 
the structure function D and the correlation matrix R a p 
are "second-order statistics" because they are explicitly 
written as second-order moments of either the full density 
V or of the reduced two-point probability density P 2 . 

In general the turbulence may not be isotropic, espe- 
cially if (B) = Bq ^ 0. It is then of interest to compute 

correlations in directions relative to Bq = Bq/Bq (and 
a perpendicular direction no). We will also consider a 
locally defined mean field in which the local mean direc- 
tion H = H/\H\ is itself a random variable that satisfies 
(H) = Bq. We have in mind several possible ways to 
determine H, for example as Hy(x) = V^ 1 J v d 3 r B, 
defined in subvolumes or "boxes" of volume V centered 
on x; or as H l = L^ 1 J L dr B, defined in co-linear sam- 
ples of length L (as used in single spacecraft estimates 
in the solar wind, under Taylor's hypothesis); or per- 
haps as H r (x) defined as the average of B at x and x' , 
\ \B(x + r) + B{x)\. These may be called, respectively, 
volume averages, line averages, and point averages. 

At this point we recall a generalization of standard 
structure functions (Milano et al. 2001) that formally 
coincides with D{r), 

D(s) = {\b(x)-b(x + sc)\ 2 ), (5) 

except that the direction c is permitted to be a vari- 
able direction, in particular a direction determined rel- 
ative to the local estimates H . For special cases, the 
separation can be chosen to be either parallel to the lo- 
cal mean magnetic field, c — » H , or else perpendicular 
to H: c — > n = e x H/\e x H\ for reference direc- 
tion e. This enables us to define the (locally) parallel 
structure function Z)H(s) = {\b(x) — b(x + sH)\ 2 ) and 

7 In the solar wind (Jokipii 1973), and wind tunnels (Monin 
& Yaglom 1971, 1975), rapid sweeping past detectors provides a 
useful connection between temporal and spatial statistics, through 
the Taylor frozen-in flow hypothesis. Here we will not be concerned 
with time statistics. 
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the local perpendicular (transverse) structure function 
D (s) = (\b(x) — b(x + sh)\ 2 ). Our main concern here 
will be to better understand the relationship between lo- 
cal and global anisotropy, and between D and D. We 
conclude that however appealing their similarity might 
be, one can anticipate fundamental differences. Notably, 
since D is defined in a stochastic coordinate system it 
is thus a conditional statistic. D is therefore necessar- 
ily a higher order statistical quantity. In contrast, D is 
defined in a fixed coordinate system and is related di- 
rectly to the second-order statistics, such as correlation 
functions and spectra. 

Milano et al. (2001) showed that structure functions 
defined as above can be employed to quantify both 
global anisotropy and local anisotropy. For separa- 
tions s lying in the inertial range and smaller, and for 
cases with nonzero Bq, global anisotropy takes the form 
D(sho)/ D(sBq) > 1, which is equivalent to enhanced 
perpendicular anisotropy seen in spectra (Shebalin et al. 
1983; Oughton et al. 1994). Local anisotropy is typi- 
cally characterized by D ± (s)/D^(s) > 1, which corre- 
sponds to stronger gradients perpendicular to a local 
preferred direction. Interestingly, the local anisotropy 
was found to be greater than the global anisotropy at a 
given scale, D- L (s)/D^(s) > D(sn ) / D(sB ) > 1 (Mi- 
lano et al. 2001). It is noteworthy that local anisotropy 
was found even in cases that are globally isotropic with 
B = 0. A related study (Cho & Vishniac 2000) found 
similar local anisotropy, and argued for a correspondence 
to anisotropic magnetic spectra in the "reduced MHD" 
regime (Goldreich & Sridhar 1995). 

2. RESULTS 

Milano et al. (2001) also suggested that there was ev- 
idence for progressively greater anisotropy as the mean 
field is computed more locally. However the available 
bandwidth of the simulations limited the strength of that 
conclusion. We present a new result here that demon- 
strates this effect with higher resolution simulations, 
employing (dealiased) incompressible 3D MHD spectral 
method simulations at a spatial resolution of 512 3 . Thus 
the maximum retained wavenumber is /c ma x = 170 (fc = 1 
corresponds to the longest allowed wavelength in the pe- 
riodic box). The runs are initial value problems, de- 
caying turbulence with initially excited wavenumbers of 
k = 1 to 6, and initial fluctuation energy normalized 
such that (b ■ b) = (v ■ v) = 1. The resistivity /j, and 
viscosity v, are selected to ensure good spatial resolu- 
tion, meaning that k max /k^i aa > 3 at all times (Wan 
et al. 2010). The dissipation wavenumber is defined as 

fcdiss = 1/Adiss = (e/^ 3 ) 1 ^ 4 , where e is the rate of (to- 
tal) energy dissipation. Initial Reynolds numbers arc 
R = l/v = Rm = 1/(1= 2000. 

We report on two runs, one with Bq = and another 
with Bq = 1. Broadband energy spectra develop and 
the data are analyzed near the time of maximum dissi- 
pation. We examined the anisotropy relative to local 
mean fields determined variously as volume, line and 
point averages. Corresponding estimates of Z)"(s) and 
D- L (s) were computed in the same way. For volume av- 
erages, or line averages with V 1 ^ 3 ~ i » A c , the cor- 
relation scale of the fluctuations, we expect that Hy 



provides increasingly refined estimates of Bq. Indeed, 
for one dimensional random process, the classical crgodic 
theorem (Monin & Yaglom 1971, 1975) guarantees that 
Hibl^Hl = (B) = B . 8 

To assess the degree of local anisotropy, the local struc- 
ture functions are computed a locally-determined mean- 
field-aligned axisymmetric coordinate system. That is, 
estimates of Z)"(s) and D ± (s) are accumulated and av- 
eraged. Several methods are used to determine the local 
mean field - volume, line and point averaging, as de- 
scribed above. Figure 1 shows the ratio D J -(s)/D^(s) 
which provides a measure of the variation of anisotropy 
with lag s, for several different determinations of the lo- 
cal mean field. Also shown is the anisotropy ratio relative 
to the global mean magnetic field. It is clear that (1) the 
degree of anisotropy is greater at smaller lags; and (2) 
that anisotropy is more pronounced when the mean mag- 
netic field is computed over a smaller region. This was 
found by Milano et al. (2001) and Cho & Vishniac (2000) 
and confirmed in numerous other studies. As has been 
reported previously, local anisotropy is present not only 
for the nonzero Bo case (Fig. lb) but also for zero global 
mean field (Fig la). We thus confirm that stronger gra- 
dients are produced perpendicular to the magnetic field, 
and also that local anisotropy is stronger than global 
anisotropy. 

Clearly the local structure functions are of great in- 
terest, and we now discuss their formal nature. It is 
convenient to introduce the generalization D(sn,s±) = 
(\b(x) — b(x + snH + s±n)\ 2 ), which makes explicit the 
possibility of arbitrary separations, both parallel (s||) to 
and perpendicular (s± ) to the local mean magnetic field. 
Evidently £>U(s) = D(s,0) and D-^s) = D(0,s). 

The quantity D can be formally obtained from the ex- 
pression for D(r) in Eq. (2) by the replacement r — > 

snH + s±h in Eq. (4). However some care is required 
in interpreting this procedure due to the new stochas- 
tic variables that appear in the arguments. (For exam- 
ple, the implication is that the differential d 3 b' becomes 
stochastic.) To avoid complications we choose to write 
the defining relation in terms of the full probability den- 
sity of realizations, that is, 



D(8 h 8 X ) 



- l r dT 



\b(x + s\\H + sxn) -b(x)\ 2 V(T). 

(6) 

This new expression is not a simple coordinate transfor- 
mation of Eq. (4) because the coordinate unit vectors 

H and h are themselves random variables, dependent 
upon r. However the new independent variables su and 
s± are non-random coordinates in the stochastically ro- 
tated reference frame of the local systems. It is clear 
that -D(s||,sj_) remains a second-order statistical quan- 
tity only if the coordinate system is fixed and H and h 
become fixed vectors, or if the probability distributions 
are insensitive to directions. The former case is obtained 
asymptotically as Bq — > oo and the mean field direction 

8 For line averages, this statistical limit is guaranteed when the 
traced correlation function tends to zero at large lags \r\ — > oo 
fast enough so that J °° dx xR(x, 0, 0) is bounded for any choice of 
Cartesian axis x. 
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Fig. 1. — Anisotropy ratio D- L (s)/D^ (s) vs. lag s at inertial 
range and larger scales. Point and line average estimates of local 
mean field H are used ("AB" averages the two points x, x' , "lag de- 
pendent" integrates from x to x' .) Perpendicular excitation is rel- 
atively enhanced, except for lags s > 1 where isotropy holds (N.B. 
energy containing eddies are initially isotropic), (a) Results from 
So = simulation, with A c = 0.34 and A^iss = 0.021. Anisotropy 
is stronger for more localized estimates of H (e.g., shorter line av- 
erages). Randomizing phases (Gaussianizing - see text), produces 
isotropic results for all averaging methods and at all lags. There- 
fore, anisotropy is associated with non-Gaussian correlations, (b) 

Results from Bq = 1 simulation, with A^ = 0.34, Ac = 0.43, and 
-^diss = 0.022. This case is anisotropic relative to the globally de- 
termined Bo i and all methods based on local determinations of the 
mean field show further enhancement of anisotropy. For this run, 
phase randomization removes the enhanced local anisotropy, but 
global anisotropy remains. 

becomes statistically sharp. The latter case is obtained 
for highly symmetric cases such as isotropic turbulence. 
In general, unless the density V is invariant under the 
stochastic changes of H, the statistical nature of D be- 
comes higher order. 

The latter point and its consequences can be made 
more explicit by considering special cases. For small 
separations purely in the parallel direction (sj_ = 
and small S||), and considering the decomposition H = 



Bq + h, we may expand 



D(s, 0) = (\sH ■ V6(as) + y-H^ V»Vj6(x) + 0(^)| 2 } 



= s 2 {\(B + h)-Vb\ 2 )+O(s 3 ). 



(7) 



Suppose that \b\/B C 1 is small enough to justify 
an additional expansion in that parameter. This yields 
D(s,0) = s 2 B 0i B 0j ((V i 6 )t )(V i 6 fc )) +O(hbb/B ). The 
first term is a second-order moment. The remaining 
terms include third, fourth, and higher-order moments. 
Only in the asymptotic limit, in which H — ^ Bo, is a 
second-order moment recovered. 



As a second special case suppose that (i) Bq = 0, 
so the global anisotropy may be zero, and (ii) that 
the local mean field estimate is completely localized so 
that H = h = b. Then once again beginning with 
the small s expansion Eq. (7), we find that D(s, 0) = 
s 2 (\b-\7b\ 2 /\b\ 2 ) + 0(s 3 ), which involves moments higher 
than second-order. Even in the very special case that 
the magnetic field fluctuation is "arc polarized" and has 
constant magnitude er, this still only reduces to a fourth- 
order quantity D(s, 0) w (s 2 /a 2 )(\b ■ V6| 2 ). 

Observing that the relationship between the structure 
function D(r) and the local field aligned structure func- 
tion D is somewhat analogous in form to the relation- 
ship between Eulerian and Lagrangian correlation func- 
tions, one might be tempted to seek similar approxima- 
tions to connect them. One possibility, based on the 
ideas underlying Corrsin's hypothesis of independence 
(see, e.g., McComb (1990)) is to treat the distribution 
of the magnetic field fluctuations as independent of the 
distribution of the mean magnetic field directions. Thus, 
symbolically one might attempt an approximation such 
as P 3 (b,6',ff) P 2 (b,b')P H (H), from which it would 
follow that £>(s||,sj_) = fdfl H C(si\H + s±n)P H (H), 

where Ph{H) is the distribution of mean field direc- 
tions and dfiij the differential of solid angles associated 
with those directions. However such an approximation 
must fail, as can be readily seen: Suppose that C(r) 
has a strong perpendicular anisotropy. Then the ran- 
dom distribution of mean field directions Pjj will dull the 
sharpness of the anisotropy by averaging parallel separa- 
tions with perpendicular separations. This independence 
hypothesis would thus produce local anisotropy that is 
weaker than the global anisotropy. This is inconsistent 
with both simulation results and solar wind results and 
therefore this approximation is invalid for the turbulence 
of interest. 

In order to demonstrate their dependence on higher 
order correlations, we examined numerically the ef- 
fect on local anisotropy of a Gaussianization or phase- 
randomization process. This procedure was carried out 
for the same simulation datasct described above. In par- 
ticular for both the Bq = and Bq = 1 cases, we modi- 
fied the Fourier coefficients by randomizing their phases 
while keeping their magnitudes unchanged. The effect 
of phase randomization is to produce a signal that is 
Gaussian, lacking coherency associated with phase cor- 
relations. However this process does not modify the en- 
ergy spectrum. Employing the phase randomized signal, 
we again compute the locally defined structure functions 
using the same set of methods for determining the lo- 
cal mean field that was described earlier. The scaling 
of the anisotropy with lag is shown in Fig. 1, where it 
is compared with the original simulation data for which 
phase coherency, if present, was maintained. The phase 
randomization has a dramatic effect in both cases: when 
there is global isotropy, the phase randomization com- 
pletely destroys the local anisotropy. When there is 
global anisotropy the phase randomization completely 
eliminates the local enhancement of anisotropy. 

3. DISCUSSION 

The issues described above appear to have immediate 
relevance to a number of analysis procedures that rely 
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on local accumulation of data. For example, even apart 
from the issue of anisotropy, the discussions presented 
here are also relevant to studies that for various reasons 
employ short datasets to define a "mean" magnetic field 
and similarly local correlation or structure functions. 9 
Invariably use of short intervals (less than A c , or its tem- 
poral equivalent) time series leads to poorly determined 
global statistics, although an interpretation in terms of 
local conditional statistics may still be meaningful (e.g., 
Sahraoui et al. 2009, 2010; Alexandrova et al. 2011). 
However, the connection to various orders of statistics 
needs to be established. 

Another popular approach is to employ wavelet tech- 
niques to determine "local spectra" as well as measures 
of local anisotropy analogous to D (e.g., Horbury et al. 
2008; Wicks ct al. 2010, 2011; Chen et al. 2011). This 
approach is also based on conditional statistics, since 
it determines a local mean field at each scale and then 
accumulates data relative to the direction of that local 
mean field. It follows that such wavelet spectra will also 
be higher (than second) order statistical quantities, and 
thus are distinct from the actual spectra. 

In conclusion, an examination of structure functions 
computed relative to locally determined mean magnetic 
field directions reveals that such quantities involve higher 
than second-order moments of the underlying probability 
distributions. This property emerges because the local 
mean field direction becomes a random variable. Con- 
sequently these structure functions are computed in a 
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stochastic coordinate system, and involve averaging over 
the magnetic field at two positions, b, and £>', and also 
the mean magnetic field direction H. Observing that 
the energy spectrum and related structure function are 
second-order moments, the fact that the locally-oriented 
structure functions involve higher order moments im- 
plies that the information they contain is not identical 
to that of the energy spectrum. The additional random 
degree of freedom docs not appear to be amenable to 
an adaptation of the Corrsin independence hypothesis, 
given that the local anisotropy is stronger than the global 
anisotropy. Gaussianization (phase randomization) of a 
turbulent field eliminates the enhancement of the local 
anisotropy, confirming its sensitivity not just to higher 
order correlations, but in particular to higher order non- 
Gaussian correlations. Evidently the phenomenon of lo- 
cally enhanced anisotropy involves fundamental physics 
that is embedded in the higher order statistics, as is the 
case for intermittency and coherent structures that are 
generated by turbulence. This information is beyond the 
scope of what can properly be described with just the 
energy spectrum. 
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9 This may be motivated by interest in local fluid physics, local 
kinetic physics, or simply lack of available high-cadence data. 



